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Abstract 

We report an experimental and theoretical study of the antiferromagnetic S=h chain subject 
to uniform and staggered fields. Using inelastic neutron scattering, we observe a novel bound- 
spinon state at high energies in the linear chain compound CuCl2 • 2((CD3)2SO). The excitation 
is explained with a mean- field theory of interacting S=^ fermions and arises from the opening 
of a gap at the Fermi surface due to confining spinon interactions. The mean-field model also 
describes the wave-vector dependence of the bound-spinon states, particularly in regions where 
effects of the discrete lattice are important. We calculate the dynamic structure factor using exact 
diagonalization of finite length chains, obtaining excellent agreement with the experiments. 



PACS numbers: 75.25,+z, 75.40. Gb, 75.10.Pq 



I. INTRODUCTION 



When quantum particles interact, they can team up to form new particles with fractional 
charge or spin as observed in two-dimensional electron gases or low- dimensional spin lattices. 
A particularly fruitful model system to study emerging new particles and their interactions 
has been the antiferromagnetic (AF) S=| Heisenberg chain. It does not order even at 
T = K due to strong quantum fluctuations and its elementary excitations are spinons 
carrying fractional S=h, which interact only weakly and are unbound particles.— Spinons 
were observed in several materials, including KCuF^, BaCu2Si 2 07^ and copper pyrazine 
dinitrate.— Strong interactions between spinons can arise from the breaking of spin rotational 
symmetry as for example in the three-dimensional coupled chain antiferromagnet where the 
mean field associated with the long-range ordered ground state restricts spin fluctuations. 
This creates an attractive potential for the spinons and at low energies they condense into 
spin-wave excitations carrying 5=1, 5 as observed in KCuFg& and BaCu2Si207.— 

Recently excitations were observed in the S— | chain antiferromagnet 
CuCl2 • 2((CD 3 ) 2 SO) (CDC) subject to uniform and staggered magnetic fields which 
were interpreted as bound spinon states.^ In the long- wave- length limit, the experiment 
demonstrated that the low-energy excitations of a S=| chain in a staggered field correspond 
to the soliton and breather excitations of the quantum sine-Gordon modeL 9 However, 
the sine-Gordon model is only valid for a very restricted range of wave-vectors and does 
not apply for excitations at smaller length scales where the discreteness of the lattice 
becomes important. More importantly, it does not fully describe the mechanism by which 
the staggered field produces an attractive potential that binds spinons into the long-lived 
dispersive S—l excitations as observed in the experiment. 

In this paper, we analyze the dispersion of the observed bound spinon states^ in detail, and 
we report the observation of a bound-spinon state at high energies in CDC. In this system 
which has a nearest neighbor exchange J = 1.5 meV, a staggered field of the order of 1 T, 
which corresponds to a Zeeman energy g/J<BH ~ 0.1 meV, qualitatively affects the excitation 
spectrum to energies more than twice J at 3.4 meV. This is in stark contrast to spinon 
binding in coupled chain magnets where the effects are only apparent for hu ~ /c^Tat.— A 
simple mean-field theory of fermions carrying 5*= | in one dimension captures the wave- vector 
dependence of the bound spinon states, and explains the high-energy excitation through 
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the opening of a gap at the Fermi surface. This model represents a first step towards a 
comprehensive description of the incommensurate excitations in AF S—^ chains subject to 
staggered fields. 



II. EXPERIMENTAL 

CDC was identified as an AF S—k chain system in which a staggered (/-tensor and/or 
Dzyaloshinskii-Moriya (DM) interaction o 10 ' 11 lead to a staggered field H st upon application 
of a uniform field H. In CDC, a uniform magnetic field H = (0, 0, H) along the c-axis 
generates a staggered field H st = (iJ st ,0,0) along the a-axis, and the Hamiltonian can be 
written as 

H = J2J$i S m - g^sHSt - gafisH^-iySf , (1) 

i 

where g c and g a are the uniform part of the gyromagnetic tensor along the c and a-axis, 
respectively. The staggered field is given by 

H st = ^DH + ^H, (2) 
g a g a 

where g s is the staggered gyromagnetic form factor and D = |D| is the length of the DM 
vector, which points along the b-axis in CDC. The nearest-neighbor spin exchange along 
the chain is accurately known from susceptibility measurements^ and inelastic neutron 
scattering.— The spin chains run along the a-axis of the orthorhombic crystal structure 
(Pnma)^ with the Cu 2+ ions separated by 0.5a ± 0.22c. Wave vector transfer is indexed 
in the corresponding reciprocal lattice Q(hkl) = ha* + kh* + lc*, and we define the wave- 
vector transfer along the chain as q = Q ■ a. Due to weak inter-chain interactions, CDC 
has long-range AF order in zero field below T/v = 0.93 K with an AF wave- vector Q m = a*. 
An applied field along the c-axis suppresses the ordered phase in a second order phase 
transition at H c = 3.9T, 10 indicating that inter-chain interactions favor correlations that are 
incompatible with the field-induced staggered magnetization.— At fields much greater than 
H c , the staggered fields thus arise mostly from a staggered (/-tensor and DM interactions 
and not from interchain interactions. 

The neutron scattering experiments were performed on 7.76 g of deuterated single- 
crystalline CDC. The measurements were carried out using the SPINS triple axis spectrom- 
eter and the DCS time of flight spectrometer at the NIST Center for Neutron Research. 



The SPINS measurements were performed with a focussing analyzer covering 7° in 20 scat- 
tering angle set to reflect Ej — 5 meV to the center of the detector. A Be filter rejected 
neutrons with energies higher than 5 meV from the detection system. The measurements 
were performed with the strongly dispersive direction along the scattered neutron direction 
to integrate over wave-vectors along a weakly- dispersive directions. The experimental con- 
figuration for the measurements made using the DCS spectrometer and the conversion of 
those data to absolute units are described in detail elsewhere.— They were performed with an 
incident energy i?j=3.03 meV and the incident beam parallel to the a-axis for configuration 
A, and with an incident energy Ei=4.64 meV and angle of 60° between the incident beam 
and the a-axis for configuration B. 

k (r.l.u.) 



-3.33 -3.08 -2.80 -2.49 -2.1 




1 2 3 4 5 

ficti (meV) 



FIG. 1: Neutron scattering intensity as a function of energy transfer for zero-field and 9.86 T 
measured using SPINS. The dashed line is a fit of two Gaussians convolved with the resolution 
function given by Cooper and Nathans^ The solid line shows the calculated intensity obtained 
from the exact diagonalization of finite chains for H = 11 T and scaled to the data, with the non- 
magnetic background given as the straight dashed line. The dashed-dotted line is the exact two- 
spinon cross-section of antiferromagnetic 5*=| chain 15 convolved with the experimental resolution 
function. Inset: Excitation energy of the higher-energy mode for H=9.86 T as a function of 
wave- vector transfer along the chain direction. 
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III. EXPERIMENTAL RESULTS 

Figure [T] demonstrates the dramatic changes that CDC undergoes upon application of 
a magnetic field. In zero field, the neutron scattering spectrum for q = OAtt consists of 
a strong peak which corresponds to the two-spinon continuum, whose band width, at this 
wave vector is narrow and barely distinguishable from the experimental resolution. In a 
#=9.86 T ~ lJ/{g(JL B ) field applied along the c-axis, which also induces a staggered field 
along the a-axis, the scattering includes two resolution-limited excitations. According to the 
quantum sine-Gordon theory, the lower-energy excitation corresponds to a bound-spinon 
state which develops into solitons and breathers at long wave-lengths. The high-energy 
excitation at 3.4 meV ~ 2.2J, however, does not have a simple interpretation in terms of 
the sine-Gordon model. Its magnetic nature is apparent on account of its field dependence. 
The inset to Fig.^shows the dispersion of this excitation, which has a maximum at g=0.47r. 

The dispersion of the high-field excitations at lower energies, Tuo < 2 meV, is illustrated 
in Fig. 121 for two different chain wave-vectors q. For q=0.7ir, there are two maxima as a 
function of energy, corresponding to the well-defined modes developing into the sine-Gordon 
soliton and breather excitations at long wave-lengths. For g=0.47r, only one peak is clearly 
observed because of the weak intensity of one of the modes in this wave-vector region. 
The dispersion of these field-induced resonant modes was determined as a function of the 
chain wave-vector, q, by fitting resolution-corrected Gaussian line-shapes to the observed 
scattering. The adjusted excitation energies are shown in Fig. Ot, illustrating that magnetic 
spectral weight generally shifts to lower energies in an applied field. However, due to the 
i7 st -induced gap, this effect is much less pronounced than in a uniform field as the ground 
state energy increases-^ 

IV. MEAN-FIELD THEORY 

We now present a simple mean-field theory of interacting S=~ fermions which captures 
both the emergence of a new excitation at high energies and the dispersion of the bound 
spinon states. As it is well known, a mean-field approach usually is not adequate to solve a 
one dimensional system. This is because, for finite range interactions, the fluctuations are 
strong enough to preclude a non-zero value of an order parameter at any fine temperature. 
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FIG. 2: Neutron scattering intensity as a function of energy transfer at 11 T for two different chain 
wave vectors, measured using the DCS instrument with an incident energy £^=4.64 meV. The 
solid line is the spectrum calculated from exact diagonalization of finite chains in absolute units, 
taking into account the polarization dependence of the experiment. 



Moreover, if the order parameter is associated with a continuous symmetry, its mean value 
is zero even at T=0. Consequently, a mean field theory that assumes a non-zero value of the 
order parameter with excitations that are originated by small fluctuations of such quantity 
cannot be a good description of a general one dimensional system. 

In the present case, the system is invariant under rotations around the z axis for H st = 
and the candidate to be the order parameter is the xy planar component of the staggered 
magnetization M±. As expected for a continuous symmetry and a gapless spectrum, the 
system is critical at T = 0, i.e., the order parameter has divergent fluctuations. However, the 
staggered field H st = couples linearly with the order parameter MJ* along the x direction. 
Consequently, for H st ^ 0, the U(l) rotational symmetry is explicitly broken and the mean 
value of Mjf becomes non-zero at any temperature. This also changes the nature of the 
excitations. The spinons (kink and antikinks) are no longer the low energy quasiparticles of 
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FIG. 3: (a) Dispersion of the excitations at H=ll T and the lower bound of the zero-field two- 
spinon cross-section^ obtained through fits to the zero-field data. The broken lines are predicted 
thresholds for continua for a 5*=^ chain in a uniform field^ The dashed-dotted line is the des- 
Cloizeaux-Pearson lower bound for excitations in S=^ chains in zero field 18 for J=1.5 meV. The 
solid lines correspond to the dispersion obtained from Gaussian fits to the spectra obtained by the 
exact diagonalization of finite chains. The dotted line corresponds to the mean-field dispersion. 
The dashed double-dot line is a guide to the eye for the high-energy mode above 3 meV, which was 
measured at 9.86 T. (b) Integrated intensity of the two resonant modes at H=ll T as a function 
of wave- vector transfer. 



the system since a pair of them is now confined by a linear potential. A similar effect occurs 
when we increase the dimension of the system due to the interchain interaction. Therefore, 
we expect a mean-field treatment to be good approach for high enough values of H st . 

For H st = 0, we know that the theory must be critical with the associated linear soft 
modes shown in FigJ3k- We also know that one dimensional spin system can be mapped into 
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a fermionic system. In addition, a non-interacting fermionic Hamiltonian has a ground state 
that is also critical and has linear soft modes like the ones shown in FigOK- Therefore, it is 
convenient to use a fermionic representation for the mean-field approach. For this purpose, 
the spin degrees of Cu 2+ ions are described in terms of fermionic creation and annihilation 
operators: 

Sj = 2 X] C }a (T aa' C ja' J (3) 
a,a' 

where v = x,y,z and a v are the Pauli matrices. Using this fermionic representation, 
Baskaran, Zou, and Anderson^ proposed a mean-field theory (MFT) to treat low dimen- 
sional Heisenberg spin S = 1/2 Hamiltonians. The MFT was generalized to SU(N) spin 
models (for large N) by Affleck and Marston.— Arovas and Auerbacb2i studied this theory 
in comparison with the Bethe Ansatz solution and showed that the fluctuation corrections 
are important in enforcing the Gutzwiller projection. Using an extended version of this 
MFT, we will study here the ground state properties and the spin dynamics of the S=^ 
chain in a staggered field given by TC. 

Apart form an irrelevant constant, the expression for TC in the fermionic representation 

is: 

W = ~ 4a C i+lA+la' C ia> ~ ° n i° ~ y^B^st ^(-1)^0 & , (4) 

with a = —a. Since there is one spin per site, there is a constraint on the fermion occupation 
number: rij = J2cr n i<r = 1 with n ia = c\ a c ia . For the Heisenberg term, we will use a 



linear combination of the mean-field decoupling introduced in Ref. 



19 and the other natural 



decoupling in the presence of a staggered field along the x direction: 



to 



(5) 

where the Lagrange multiplier or chemical potential, A, enforces the constraint of one spin per 
site at mean- field level. We are assuming translational invariance for 7^ 7 = Z)<r( c L c j+icr) = 
7i, and a staggered dependence for the effective field Si = J2a( c ia c ia) — S(— 1)\ This 
staggered dependence is induced by the field H st . In momentum space, this leads to 



Hmf = ( ~Jl cos(£;) - ^-g c fj, B H) 4 a c ka - \{g a ^BH st + J5)(4 +WCT c ka + c^c k+Tl 



-7r<fc<7T,cr 



(6) 
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which can be written in the matrix formulation as 



H M F{k) a 



-Jjcos(k) - a\g c ^ B H -\{g a ^BH st + J 8) 
-\{g a ^BH st + J 5) J7 cos(fc) + a\g c [i B H 

for — 7r/2 < k < 7r/2. The eigenvalues of this matrix, 



<~to 71 ±\/ ( Jjcos{k) + a^fi B H\ + ^(g a ^BH st + J8) 2 



are the energies of the quasi-particle operators, 



with 



Vka 



t\ a + J7 cos(fc) + a\g c {i B H 



e\ c + J7 cos(fc) + a\g c \x B H] + \{g a ii B H st + J 8) 



\g a ^BH st + J8 



[e~£ a + J7 cos(&) + o\g c ii B H} + \(g a fi B H st + J8f 
Hmf is diagonal in the new basis: 



HMF = (4*Pk*Pk* + e k* a k* a ka) ■ 

-7r/2<fc<7r/2,cr 

The mean-field parameters 7 and 8 are given by the self-consistent equations: 
1 r 71 "/ 2 

7 = 77- E / , COS ( fc ) ( M L - U L) [(«L«fca) - (PkJka)} d k , 
2,71 J-tt/2 



1 r 71 "/ 2 

All „ J-tt/2 



(7) 



(8) 



(9) 



(10) 



(11) 



(12) 



The value of A is determined by imposing the average occupation per site to be equal to 1: 

rir/2 



r' K«L«jJ + (PlJka)]dk = 1 

ATI p J-vr/2 



(13) 



When H = H st = 0, the integration over the phase fluctuations of the local field 7^ 
renormalizes the value of 7 given by Eq. (fT2"j) : 7 = This renormalization improves 

considerably the comparison with the exact des-Cloizeaux- Pearson^ two-spinon threshold 
{lex — 7r /2). To improve the quantitative comparison of the MFT with the experiment 



and the exact diagonalization results, we will assume here that the same renormalization 
factor, 7r/v^2, must be applied when H and H st are finite. The value of 7 is a measure 
of the effective strength of the spin fluctuations introduced by the Heisenberg term. More 
specifically, 7J is the Fermi velocity that in the original spin language corresponds to the 
velocity of the spinon excitations. 

Fig. |3] shows the evolution of the fermionic bands when the uniform and the staggered 
fields are applied. In the absence of these magnetic fields (Fig.Hk,), there is only one band and 
the two-spinon cross section is associated with the continuum of particle-hole excitations. 
The dispersion relation for the lower branch is 7JI sin(g)|. When a uniform magnetic field 
H 7^ is applied (Fig. 0Jd), the spin up and down bands are split by the Zeeman term. As a 
consequence, there is a change Sq of Fermi wave vectors \qp\ = ir/2±Sq and a corresponding 
change in the wave vectors of the zero energy modes: the energy of the transverse modes goes 
to zero at q = 7r and q = 25q, while the longitudinal excitations have gapless modes at q = 
and q = re — 25q. From Eq. (JBJ), the main effect of a non-zero staggered field H st is to open 
a gap at the Fermi level, i.e., the fermionic system becomes an insulator and the spectrum 
is gaped for any excitation (Fig. Ht). The gap results from the inter-band scattering which 
is introduced by the staggered field H st . According to Eqs. Q, the degree of mixing is 
maximum at q — qp. The emergence of this gap is consistent with the experimental data 
shown in Fig. El 

To study the excitations of the new ground state induced by H st it is necessary to calculate 
the neutron scattering cross section within our MFT. The neutron scattering cross section 
for the transverse excitations {y — x,y) at T = K is given by: 



where 77 = — 1 for v — x, T) — 1 for v = y and — it < q < 71. Note that in Eq. (jl4J) . k + q 
must be contracted to the reduced Brillouin zone (k + q = k + q + mr). The cross section 
for longitudinal excitations is: 




(14) 



S ZZ (q, u) = — J2i / {u k+ qaU ka + V k+qaVk(T ) 2 5(uJ - 4 +q ^ + £far ) 




(15) 
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FIG. 4: Schematic representation of the Fermi particle dispersion in (a) zero field, (b) uniform field 
H and (c) uniform field H and staggered field H st with H > H s t. In zero field the spin up and 
down particles with wave-vector k have the same dispersion. 

Equations (fT4"j) and (fT5j) reveal well-defined excitations that are determined by the can- 
celation of doo/dk, i.e., the divergence of the Jacobian. In Fig. EH, we show the different 
branches of the transverse excitations. The lower branch corresponds to the dispersion rela- 
tion of the low energy excitations. In agreement with the experiment (see Fig. EJ), there are 
two minima, one is located at the incommensurate wave vector qj = 2 aicsm(g c HBH/2~/J) 
and the other one occurs at q = tt. It is interesting to note that for v = x the intensity 
of this branch goes to zero at qi = 25q due to a cancelation of the matrix element that 
multiplies the delta function in the integrand of Eq. (fl4j) . The black and the blue curves 
of Fig. EH are the upper boundaries of inter-band particle-hole excitations associated with 
the transverse modes (Fig. 0Jd). More specifically, the black curve results from excitations 
in which an electron is annihilated in the lower band and created in the upper band, while 
for the blue curve the process is the opposite. The green curve of Fig. Eb is the upper 
boundary for the intra-band particle-hole excitations that describes the longitudinal modes. 
Note that these boundaries appear with dashed lines in the spectrum of excitations with the 
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other polarization (Fig. HJ. This is a consequence of the inter-band q = it scattering which 
is introduced by the staggered field H st . The dashed line just indicates that these "shadow" 
branches have a very small intensity. 




q/71 q/71 

FIG. 5: a) Transverse and b) longitudinal excitations obtained from the mean- field theory, as 
described in the text in more detail. 

The gap in the low energy spectrum is not the only qualitative change introduced by 
the staggered field H st within the mean-field approach. Since the Fermi wave vectors qp = 
±ii/2±5q are now extremal points of the new bands (see Fig. Eh), we expect the emergence 
of a new branch of transverse excitations associated with transitions between points that 
are close to qF = z F'k/2 ± Sq (qF = ±7r/2 ± 5q) and points in the proximity of q = 
(q = tt). This new branch of excitations (see the red curve in Fig. ®l) coincides with the new 
excitation which is experimentally observed at high energies (see Fig. and thus explains 
our experimental results. For the longitudinal polarization, the new branch of excitations is 
shifted by it relative to the transverse polarization (see the red curve in Fig. EJd). Since the 
maximum of this longitudinal branch is located at q = it/ 2 + 5q, it is difficult to distinguish 
this branch from the breathers in the experimental data. The energy of this maximum 
is close to 3meV according to the MFT while the experimental value is 3.4meV. This is 
reasonable given that the MFT is not an adequate approximation to give a quantitative 
description of the excitations. 

V. EXACT DIAGONALIZATION OF FINITE LENGTH CHAINS 

The intensities of the different branches are not properly described by the mean-field 
equations (fTl| and For instance, the MFT predicts a high intensity for the upper 
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boundary of the two-spinon excitations even at zero field for which very accurate calculations 
are available. 16 This is clearly an artifact of the MFT. To obtain an accurate description of 
the intensities and the energies of the different branches we complemented our analytical 
approach with the exact diagonalization of finite size chains. Using the Lanczos method, we 
obtained the exact ground state of H for finite chains of length L = 12, 14, 16, 18, 22, 20, 24. 
Having the ground state, we computed the dynamical magnetic susceptibility, x{ u >q)i f° r 
all the possible wave vectors q = 0, 2rr/L, 2n(L — 1)/L of a chain of length L using the 
method introduced in Ref. 22. The wave vector , = * is present in all of the considered 
chains. The small changes in the calculated xi^i ^) as a function of L indicate that the finite 
size effects are small for the considered problem. In general, smaller finite size effects are 
expected for systems that have an excitation gap because the spin-spin correlation length is 
finite. 

The T = K structure factors were calculated using 

S aa (q,uj) = -x" aa (q,uj). (16) 

7T 

The energy spectra were obtained by convoluting the discrete spectra for finite chains with 
Lorentzian functions with a full width at half maximum, 2T = O.lmeV, in order to model 
the experimental energy resolution. The intensity of the calculated structure factors is given 
for a chain of L spins and normalized so that J2q,a I duS aa (q, uj) = S(S + 1) as required by 
the total scattering sum rule. 

The calculated neutron scattering for all chain lengths was averaged and the dynamic 
structure factor for the three different polarizations is shown in Fig. [Hjas a function of wave- 
vector transfer and energy transfer for H st =0.075H on an absolute scale. The structure factor 
S zz polarized along the uniform field contains a well-defined excitation with a minimum gap 
energy at an incommensurate wave-vector. The structure factors S xx and S yy polarized 
perpendicular to the uniform field contain well-defined excitations whose dispersion has a 
minimum at the antiferro magnetic point, with S yy having an excitation at lower energy 
than S xx . These excitations correspond to the first and second breather of the quantum 
sine- Gordon model. 

Fig. IH1 also provides evidence that the excitation spectrum contains a substantial amount 
of continuum states, as observed in the experiment £ These states lie higher in energy than 
the well-defined low-energy excitations and extend to high energies. At H st =0.075H our 
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FIG. 6: The dynamic structure factors S aa (q,u>) obtained through exact diagonalization of finite 
chains for H s t=0.075H , where a=x, y and z are the spin polarizations. The uniform field is along 
the c-axis (z polarization), the staggered field along the a-axis (x polarization). 

numerical calculations yield an expectation value of the Hamiltonian per S = | of TC = 
— 0.34J, increased from the ground state energy at zero field, H = (| — In 2) J = —0.44 J, by 
0.1J. 

The numerical data were binned and Gaussian fits were used to obtain excitation energies 
as a function of wave- vector. These results are shown in Fig. as a solid line, showing that 
the numerical calculations reproduce the dispersion relation of the low energy modes. We 
calculated the dynamic structure factor for the DCS experiment taking into account wave- 
vector dependent mixing of the polarized dynamic structure factor S aa . Fig. [7| directly 
compares the calculated and measured intensities on an absolute scale, showing that there 
is excellent agreement between the numerical calculations and the experiment. 

In addition, Figure El shows that a new branch of transverse magnetic excitations with 
the maximum around q = OAtc and Tiuo = 3.5meV emerges when the staggered field is 
present. This provides a quantitative explanation for the high energy peak that appears in 
the neutron scattering data (see Fig.^). Muller et a/.— showed that this branch is not present 
for H st = by using macroscopic selection rules. For non-zero H st , the total spin S and its 
projection along the z-direction, S z , are not good quantum numbers anymore, allowing the 
emergence of this new excitation. Our mean-field approach is in good agreement with this 
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result. 



VI. CONCLUSIONS 



In summary, we have performed neutron scattering experiments, numerical calculations 
and an analytical study investigating the AF S=^ chain in uniform and staggered fields. We 
found that the incommensurate bound-spinon states are well described by a mapping to an 
interacting fermionic model after a renormalization of the energy. The model also explains 
the emergence of a new excitations with the application of a staggered field. Our results 
suggest that the proposed mapping is more powerful than initial results suggested, and that 
it may also be useful for other quantum spin systems with a relatively short correlation 
length. 
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